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Abstract: We define sets of orthogonal polynomials satisfying the additional constraint of a 
vanishing average. These are of interest, for example, for the study of the Hohenberg-Kohn functional 
for electronic or nucleonic densities and for the study of density fluctuations in centrifuges. We 
give explicit properties of such polynomial sets, generalizing Laguerre and Legendre polynomials. 
The nature of the dimension 1 subspace completing such sets is described. A numerical example 
illustrates the use of such polynomials. 


I. INTRODUCTION 

Generalizations of Hermite polynomials H n were recently [1] proposed to describe, for instance, density pertur¬ 
bations constrained by a condition of matter conservation. Because of the constraint, such polynomials cannot form 
a complete set, but span a subspace well suited to specific applications. In particular, the polynomials T ra used in [1] 
were motivated by the consideration in nuclear physics of the Hohenberg-Kohn functional [2] and similar functionals 
along the Thomas-Fermi method [3,4]. Indeed, in such approaches, the ground state of a quantum system is shown to 
be a functional of its density p{r), and there is a special connection between p{r) and the mean held u{r) driving the 
system. It was thus convenient to expand variations of p in a basis (w m (r)} of particle number conserving components, 
5p(r) = Sp m w m (r), with the term-by-term constraint, Vm, f drw m (r) = 0. This spares, in the formalism, the 
often cumbersome use of a Lagrange multiplier. Simultaneously, it was convenient to expand variations of u in a basis 
orthogonal to the hat potential, because, trivially, a hat Su , as just a change in energy reference, cannot influence the 
density. The same basis can thus be used for 5u{r) = <5it n w n (r), since the very same condition, f drw n (r) = 0, 
induces orthogonality to a constant Su. Because of the nuclear physics context of [1], harmonic oscillators shell models 
were considered and the basis contained a Gaussian factor, e~^ r . 

The same functional approaches [2-4] are also of a general use in atomic and molecular physics, where Gaussian 
weights would be clumsy and radial properties are best htted with simple exponential weights [5]. Furthermore, in 
[1], the discussion was restricted to one dimensional problems. In the present note, we want to include two and three 
dimensional situations. We shall thus use weights of the form e ~^ r , with 0 < r < oo, but integrals will carry a factor 
r v , with v a positive exponent, suitable for dimension d. This will lead to generalizations of Laguerre polynomials. 

This note is also concerned with compact domains, of the form 0 < r < 1 for instance. This might correspond 
for instance to expansions of density fluctuations in cylindrical vessels used for chemical processes, where mass 
conservation is also in order, or maybe in centrifuges. Radial integrals with factors r and r 2 in both the constraint 
and orthogonalization conditions will lead to generalizations of Legendre polynomials. 

For any positive weight p(r), and any dimension d , a constraint of vanishing average, f drr v p,(jr)r n {r) = 0, is 
incompatible with a polynomial T of order n = 0. Therefore, in the following, the order hierarchy for the constrained 
polynomials runs from n = 1 to oo, while that for the traditional polynomials runs from 0 to oo. We study in some 
generality the “Laguerre” case in Section II. In turn, the “Legendre” case is the subject of Section III. A brief Section 
IV discusses possible applications to the study of density fluctuations in centrifuges. Section V answers a question 
which was omitted in [1], that of the nature of the projector onto the subspace spanned by the constrained states and 
the nature of the codimension of this subspace. A numerical application is provided in Section VI. A discussion and 
conclusion make Section VII. 


II. MODIFICATION OF LAGUERRE POLYNOMIALS BY A CONSTRAINT OF ZERO AVERAGE 


In this Section we consider basis states carrying a weight e~? r , in the form w n (r) = e~? r G d (r), where G^ is a 
polynomial. It is clear that Gq cannot be a finite, non vanishing constant if the constraint, / 0 °° dr r d_1 e~ 5 r Gg(r) = 0, 
must be implemented. Hence set integer labels m > 1 and n > 1 and define polynomials G d by the conditions, 


dr 


d-i e -r 


G d m (r)G d n (v)=g d 6 ri 


dr r d 1 e '• 


G d n (r)= 0, 


(1) 


1 


where S mn is the usual Kronecker symbol and the positive numbers g d are normalizations, to be dehned later. 

It is elementary to generate such polynomials numerically, in two steps by brute force, namely i) first create “trivial 
seeds” of the form, s d (r) = r n — ( r n )d , where the subtraction of the average, (r n )d = 2 n (d — 1 + n)\/(d — 1)!, ensures 
that each trivial seed fulfills the constraint, then ii) orthogonalize such seeds by a Gram-Schmidt algorithm. The first 


polynomials read, 






G\=r- 2, 

G\ = r 2 - 5r + 2, 

G\ = ? 

- 10r 2 + 20r - 8, 

G\ = r A - 17r 3 + 78r 2 - 108r + 24, 

(2a) 

G? = r - 4, 

G\ = r 2 - 8r + 8, 

G\ = r 3 

- 14r 2 + 44r - 32, 

Gl = r 4 - 22r 3 + 138r 2 - 288r + 144, 

(2b) 

G\ = r- 6, 

G\ = r 2 - llr+ 18, 

CO 

II 

CO CO 

C5 

- 18r 2 + 78r - 84, 

G 3 4 = r 4 - 27r 3 + 216r 2 - 606r + 468. 

(2c) 


All these are defined to be “monic”, namely the coefficient of r n is always 1. For an illustration we show in Figure 1 
the new polynomials G\ and Gf, together with Laguerre polynomial L\. The same Fig. 1 also shows GG 2 and L 2 . 
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FIG. 1. Comparison of Laguerre polynomials L 1 , L 2 (full lines) with new polynomials G \, G\ (long dashes), G 2 , G| (dashes). 


Rather using the Gram-Schmidt method, we find it easier, and more elegant, to generate the polynomials G 
starting from the initial table, Eqs. (2a,2b,2c), by means of the following recursion formula, 

G d n (r) = (r — d) G^r) -2rG dl _ 1 (r) + (n + d - 1) (n - 2) G d n _ 2 (r), (3) 

where the prime denotes the derivative with respect to r. Its simple structure can be proven analytically as follows: 

i) Let us first create some kind of a “less trivial seed” at order n, assuming the polynomial G^_ x is known. For 
this, try rG d _i- By partial integration, we see that, 


/»oo /*oo 

/ dr r d ~ l e - 5 r [rG'^_ 1 (r)] = 2 / dr e~i r [r rf G^_ 1 (r)] / , 

Jo Jo 


where again a prime means derivation with respect to r. Thus a d = ( rG d _ t — 2 rG d '_ 1 — 2dG^_ 1 ) makes indeed a 
less trivial seed, compatible with the constraint. Notice that the order n of this seed polynomial a d comes from the 
term rG c j l _ 1 only, the other two terms having order n — 1. Notice again that, in the table, Eqs. (2), all polynomials 
G d are monic. We can define G d as monic, systematically. Since the product rG d _ 1 respect this “monicity”, and 
since a d fulfills the constraint, we conclude that cr d is a linear combination of G d , with coefficient 1, and of all the 
lower order polynomials G^,, with 1 < m < n, but with yet unknown coefficients. 
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ii) It turns out that such coefficients vanish if m < n — 2. Indeed, an integration of a d against Gj^, weighted by 
r d ~ 1 e~ r , gives, by partial integration of the G d '_ 1 term, 

Jo°° dr e~ r r d ~ l cr^(r) G d m (r) = / 0 °° dr r d - 1 e" r [(r - 2d) G^-^r) - 2rG d, _ 1 (r)\ G^(r) = 
f 0 x dr e" r r d - 1 G^_ 1 (r) (r - 2d) G^(r) + 2/ 0 °° dr G^r) [e- r r d G d m {r )\' = 

Jo°° dr e~ r r d ~ l G d _ 1 (r) [~a^ +1 {r) - 2dG^(r)\ . (5) 

In the bracket [ ] in the last right-hand side of Eq. (5) the seed cr^ +1 has order m + 1 and, by definition, G^ is 
of order m. By definition also, G d _ ± , of order n — 1, is orthogonal to all those polynomials of lower order, that are 
compatible with the constraint. This integral, Eq. (5), thus vanishes as long as m + 1 < n — 1. It can be concluded 
that the difference, a d — G^, contains only two contributions, namely those from G d _ 2 and G^ l _ l . Explicit forms for 
their coefficients are obtained by elementary manipulations, leading to Eq. (3). Elementary manipulations also give, 

2 r G d " - (r - 2d) G dl + n G d = (n - 1) (n + d) G d _ 1 . (6) 

Here, in the same way as a prime means first derivative with respect to r, we used double primes for second derivatives. 
Finally the normalization of the polynomials is obtained easily as, 

dr e~ r r d_1 [G d (r)} 2 = (n — 1)! (n + d )\. (7) 



III. MODIFICATION OF LEGENDRE POLYNOMIALS BY A CONSTRAINT OF ZERO AVERAGE 

Legendre polynomials, and their associates and generalizations (Gegenbauer, Chebyshev, Jacobi) are defined with 
respect to the [—1,1] segment. Exceptionally in the literature, one finds shifted Legendre polynomials, adjusted to the 
[0,1] segment. We are here interested in applications to radial densities in cylinders, or the small circles of toruses, 
or spheres. Hence we shall use 0 < r < 1. It is clear that the case, d = 1, does not make an original problem, since 
Legendre polynomials, whether translated and/or scaled or not, already average to 0 as soon as their order n is > 1. 
We keep the case, d = 1, for the sake only of completeness and in this Section we consider d = 1, 2,3, with a geometry 
factor r d_1 . The weight is ^{r) = 1, hence our states are described by just a polynomial Q d of order n. It is again 
obvious that Q d cannot be a non vanishing constant if the constraint, dr r d ~ x Q d (r) = 0, is implemented. Hence 
set m > 1, n > 1, and define polynomials Q d from conditions, 

C drr d - 1 Q d m (r)g d {r)= 1 d Jmn, [ dr r*" 1 G d (r) = 0, (8) 

Jo Jo 

where the normalizations j d are again to be defined later. It is obvious that the shifted (and shrunk) Legendre 
polynomials C n (2r — 1), n > 1, satisfy both constraint and orthogonality relations for d = 1, because they are 
orthogonal to any constant polynomial, of order 0. The polynomials = C n (2r — 1) thus make nothing new. We 
turn therefore to d = 2 and d = 3, with a brute force construction as in the previous Section. But the defining 
conditions, Eqs. (8), show a difference with Eqs. (1): both orthogonality and constraint conditions now carry the 
same weight, namely fi 2 = n, while in the previous case, Eqs. (1), there were different weights, because of the 
exponentials e~ r and e~^ r . A similar difference between /i 2 and fi happened in the “Hermite” case, naturally. 

Hence now, in this Legendre case, we can Gram-Schmidt orthogonalize even more trivial seeds r", without subtrac¬ 
tions, and accept those orthogonal polynomials with order m > 1. The table of first results reads, 

Q\ = 2r — 1, Q\ = 6r 2 -6r + l, Q\ = (2r - l)(10r 2 - lOr + 1), Q\ = 70r 4 - 140r 3 + 90r 2 - 20r + 1, (9a) 

Q\ = 3 r - 2, G\ = 10r 2 - 12r + 3, Q\ = 35r 3 - 60r 2 + 30r - 4, Q\ = 126r 4 - 280r 3 + 210r 2 - 60r + 5, (9b) 

gf = 4r — 3, g\ = 15r 2 - 20r + 6, g\ = 56r 3 - 105r 2 + 60r - 10, = 210r 4 - 504r 3 + 420r 2 - 140r + 15. (9c) 

Easy, but slightly tedious manipulations validate the following recursion relations, 

n G\ = (2 n - 1) (2r - 1) 5 4 _i - (n - 1) 0 4 _ 2 , 

(n + 1) (2 n - 1) g 2 n = 2 [ (4n 2 - 1 )r - 2n 2 ] g 2 _ ± - (n - 1) (2 n + 1) g 2 n _ 2 , 

n 2 (n + 2) t/ 3 = (2 n + 1) [2n(n + 1 )r - ( n 2 + n + 1) ] g^-i - (n - 1) (n + l) 2 t/ 3 _ 2 . 
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(10a) 

(10b) 

(10c) 



and the differential equation, 


r{r-l)Q*" + [(d+l)r-d\Qi' -n[n + d)Qi = Q. (11) 

Finally the normalization of the polynomials reads, 

pOO 

7n = / dv r d_1 [Gn(r)] 2 = 1/ (2n + d). (12) 

JO 

We show in Figure 2 the plots of for n = 1, 2 and rf = 1,2,3. 


CstLeg 



FIG. 2. Modified Legendre polynomials <Ji, (full lines), , Q\ (long dashes), Qf, Q\ (dashes). 


IV. POLYNOMIALS FOR CENTRIFUGES 


The case of centrifuges is worth a short comment. As soon as the matter under centrifugation is compressible, the 
density becomes much larger at the outer edge, r = 1, than at the rotation axis, r = 0. Let h be the height of the 
centrifuge. Assume, for the sake for the argument, that one studies fluctuations about a reference density of the form, 
P{r) = Pee 


Kr z 


where the parameter K contains all informations about the angular velocity, compressibility, etc. of 

1 1 — i —i 

' ' N = M h~ 1 2 K \e K — l] , ensures the conservation of the mass 


h f 1 dr rp(r) 


the process. The factor, p c = M 
M included in the vessel. If a cause for fluctuations of p is an instability of K, the first order for density change is, 


dp ^ r,Kr 2 \e K — l] + e K — Ke K — 1 Kr 2 
dK [r) ' ' P - if 


d P t\ n 
drr -Q^( r ) = 0 , 


to 


(13) 


namely a polynomial of order 2 multiplied by e Kr . Higher derivatives with respect to K will generate similar, even 
order polynomials, with the same property, f^ drrd n p/dK n (r) = 0. An orthogonalization, under a metric oc e 2Kr , 
might be useful. This new set of polynomials will depend on K , however, since r is already scaled to a radius 1 for 
the cylinder and thus K cannot be scaled away. Because of this K dependence we do not elaborate further on this 
issue. For a large list of ad hoc polynomials and integration weights, see [6]. 
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V. PROJECTOR ON THE CONSTRAINED SUBSPACE 


For the sake of the discussion and short notations, set first d = 1, fi(r) = e 2 r , and temporarily include normaliza¬ 
tion to unity factors into both the Laguerre polynomials L n and the constrained G\. This summarizes as, 


/»oo /»oo /»oo 

/ dr [n(r)} 2 L m (r) L n (r) = Smn, / dr [n{r)f G^(r) G*(r) = S mn , / dr/x(r) G*(r) = 0, 
Jo Jo Jo 


(14) 


Then the kets and bras defined by (r\w n ) = (w n \r) = w n (r) = fi{r) G„(r) and (r\z n ) = (z n \r) = z n (r) = /i(r) L n (r) 
provide two “truncation” projectors, Vn = Y^n=i l u ’™}( u, ™| and Qn = o l 2 n)(^n|i available for subspaces where 
polynomial orders do not exceed TV. Their respective ranks TV and TV + 1, and the embedding and commutation 
relation, [Vn, Qn\ = Vn, are obvious. Obvious also is the limit, limjv^oo Qn = 1. The role of the rank one | cr/v) (er n\ 
difference Vn — Qn is to subtract from any test state, |r) = Yln=o T n\ z n), that part which violates the condition of 
vanishing average. We shall show that the elementary ansatz, 


JV 


'2 N 


I Tv) = 


\rn=0 


poo 

y^iZn) \z n ), {Zn)= dr{r\z n ), 

n—0 


defines the proper “subtractor” operator |a/v)(cr/v|. Indeed, from 


N 


N 


-1 


( Qn ~ kiv>(^iv|) | r) = | z n ) - 


n—0 


\m =0 


' N 

E< 

, n—0 


one obtains 


( N 

E( 

\ p =o 


'Pi 'p i > 


f oo N / AT \ - 1 / JV 

/ dr (r| (Qn — |crjv)(c r JV|) |t) = r n (z n ) — I (z m } 2 J ( E^") ^ 

^ n=0 \m =0 / V n =0 


N 


E^p) t p ) = °- 

Vp=o / 


(15) 


(16) 


(17) 


Hence Qn — Io'jv) (ctjv | is the projector Vn- Incidentally, the Laguerre result for cjn is very simple, because (z m ) = 
2, Vm. But the ansatz for <jn, Eq.(15), generalizes to all cases. For instance, with Hermite polynomials, odd orders 
already satisfy the constraint when integrated from — oo to oo, naturally, and thus do not contribute to <jn- Even 
orders contribute, and it is easy to verify, upon integrating from —oo to oo again, that ( Z 2 P ) 2 = n 2 2 1 ~ p (2 p — l)!!/p!. 

It may be pointed out that the condition, f dr n(r) f(r) = 0, for functions / orthogonalized, like our polynomials, by 
a metric [/r(r)] 2 , might be interpreted as an orthogonality condition, J dr f(r) [/r(r)] 2 g(r) = 0, with g(r) = [/x(r)] _1 . 
This makes g a candidate for the subtractor form factor a. This is of some interest for the centrifuge case, where a 
state function such as, for instance, e~ Kr , remains finite when 0 < r < 1. But there is little need to stress that, when 
the support of n extends to 00 , then /z -1 does not belong to the Hilbert space and cannot be used for cr. 

More interesting is the limiting process, N —» 00 , as illustrated by Figures 2-5. Figs. 2 and 3 show the shapes, in 
terms of r, of (2|'PAr|r) and (10|'P/v|r), respectively, when the projectors are made of the modified Laguerre polynomials 
G^- The build up of an approximate (5-function when N increases from N = 50 (short dashes) to N = 100 (long 
ones) and N = 150 (full lines) is transparent, although the convergence is faster when peaks are closer to the origin, 
compare Figs. 2 and 3. The slower convergence in Fig. 3 is due to the cut-off imposed by exponential weights as 
long as N is finite. Given TV, there is a “box effect”, the range of the box being of order ~ TV. A similar build up is 
observed for our other families of constrained polynomials, with slightly different details of minor importance such as, 
for instance, a box range ~ v/TV for the Hermite case. 

The box effect is even more transparent in Figs. 4 and 5, which show the shapes of subtractors (10|crjv)(o'jvk) and 
(0|erjv)(<7jv|r) deduced from constrained polynomials of the Laguerre (Fig. 4) and Hermite (Fig. 5) type, respectively. 
(For graphical convenience, the polynomials T„ and H n used for the Hermite case, Fig. 5, are tuned to a weight e~ r 
rather than e~i r , but this detail is not critical.) 

It seems safe to predict that, given an effective length A(TV) for the box, the wiggles of the subtractor will smooth 
out when TV —» 00 and that only a background ~ — 1/A(TV) will then remain. 
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pro j2 



FIG. 3. Shapes of projectors made of polynomials G„. Full line, (2|'Pi5o|r}, long dashes, (2|'Pioo|r), short dashes (2|Pso|r). 


projlO 



FIG. 4. Shapes of projectors made of polynomials G„. Full line, (10|'Pi5o|r), long dashes, (10|'Pioo|f), short dashes (10|'P5o|r). 
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FIG. 5. Subtractors made of G Shapes centered at r = 10. Short dashes, N = 10, long dashes, N = 20, full line, N = 30. 



FIG. 6. Subtractors made of F„. Shapes centered at r = 0. Stronger wiggles, shorter cut-off, dashed line, N = 50. Weaker 
wiggles, larger cut-off, full line, N = 100. 


















VI. ILLUSTRATIVE EXAMPLE: TRAJECTORIES IN DENSITY SPACE 


We return here to the toy model discussed in [1] and the corresponding, modified Hermite polynomials. The model 
consists of Z non interacting, spinless fermions, driven by a one dimensional harmonic oscillator Hq = |(— d 2 /dr 2 +r 2 ). 
The ground state density from the Z lowest orbitals reads, p(r) = Let i = 1, ..,Z and I = Z + 1, ...oo 

label “hole” and “particle” orbitals, respectively. Add a perturbation Su(r) to the initial potential r 2 /2. The first 
order variation of the density is, 


5p{r) =2^ V’iWV’/W 

Z ' l—l 


il 


If we expand 6u and Sp in that basis {w n } provided by the new polynomials, the formula, Eq. (18), becomes, 

5p m = 2 E U J—J ^ >n 11 ^ Un ’ u ~ j Wn ^ ( r ) ’ 


(18) 


(19) 


il n 


where V denotes both a particle-hole matrix element of a potential perturbation and the projection of a particle- 
hole product of orbitals upon the basis {w n }. In [1] we briefly studied the eigenvalues and eigenvectors of this 
symmetric matrix, T = T> [Eq — ifo) -1 where (Eq — Hq )^ 1 is a short notation to account for the denominators 
and the particle-hole summation, and the tilde indicates transposition. It is clear that the invertible T represents 
the functional derivative 8p m /5u n and is suited for infinitesimal perturbations. We shall now take advantage of the 
representation provided by {w n } to study finite trajectories p(u). 

For this, we consider a variable Hamiltonian, Tt m {X) = Hq + A w m (r), made of the initial harmonic oscillator, but 
with a finite perturbation A u along one “mode” w m . It is trivial to diagonalize Tt m { A) with an excellent numerical 
accuracy and thus obtain, given Z, the ground state density p(r, A). Then it is trivial to expand the finite variation, 
A p = p(r,\) — p(r,0), in the basis {u>„}. This defines coordinates Ap n (X;m) for trajectories, parametrized by the 
intensity of the chosen mode in for Au. 


rho_n 


0.1 



- 0.2 



u4 


FIG. 7. Coordinates of the perturbation density A p created by a perturbing potential Au = A 4 W 4 . Full line: 2Ap2- Long 
dashes: Ap 4 . Moderate dashes: 2 Ape. Short dashes: 4Aps- Very short dashes: 8 Apio. 



rho_n 



FIG. 8. Same as Fig. 7, but now Ait = XeWe . Full line: 4 Ap 2 - Long dashes: 2 Ap 4 . Moderate ones: Ap&. Short ones: 2 Aps. 
Very short dashes: 4 Apio. 



FIG. 9. 3D trajectory in density space. Ap 4 , Ap§ and Apg taken from Fig. 7, the latter two coordinates blown 4 times. 

In Figures 7 and 8 we show, with Z — 4, results from TLa = Hq + A 4 2 (27t)“3 15“^ ( 8 r 4 — 14r 2 + 1) e~ r and Ti.6 = 
H 0 + X e (2n)~i 105~^ (32r 6 — 128r 4 + 94r 2 —11) e~ r2 , respectively. The case, 7 i, 2 = H 0 + X 2 2 ( 2 / 7 r)i 3~3 (2r 2 —1) e -1 " 2 , 
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makes almost a harmonic oscillator and is probably of academic interest only; anyhow we verified that its confirms 
the results with H 4 and He- We use a basis {u>„} containing a factor e~ r rather than e~^ r to better match the 
same factor e~ r created by products of harmonic oscillator orbitals in the calculation of matrix elements (z p \A'u\z q ), 
but this technicality is not important for the physics. 

The main result to be observed seems to be the lack of “collectivity” for such modes and for such elementary 
Hamiltonians. Indeed, for A 4 = 2, the first five coordinates of A p read {0.016,-0.267,-0.055,0.023,0.018}, with a 
strong dominance of A/34, while for Xe = 2, these read {—0.013, —0.041, —0.376,0.008,0.040}, with a strong dominance 
of Ape. To clarify Figs. 7 and 8 , we had indeed to blow up each A p n by a factor where m is the index of the 

driver mode in potential space. Other modes than m = 4 and m = 6 show the same property: in the density space, 
a trajectory driven by A u = A w m stays close to the same w m axis in that density space, although curvatures effects, 
while somewhat modest, are not absent. Such non linearity, slight curvatures are seen in Figs. 7-8, and also in Figure 
9, where the three A/ 34 , Ape , Asets of data shown by Fig. 7 are converted into a parametric plot for a trajectory. 
For graphical purposes again, Ape and Aps are blown up 4 times to create Fig. 9. It can be concluded, temporarily, 
that the “flexibility” matrix T is not too far from being diagonal in the {rc n } basis, or in other words, that the w n 
modes indicate an approximately natural hierarchy in both the potential and the density spaces. 

A subsidiary question pops up: that of the positivity of p. Indeed, while the space of potentials is basically a linear 
space, with arbitrary signs for u(r) when the position r changes, densities p(r) must remain positive for every r. This 
creates severe constraints for any linear parametrization of A p in terms of the basis {»„}. In our toy model, it turns 
out that p(r, 0) = ( 8 r 6 — 12r 4 + 18r + 9) e~ r / 6 . Hence, if we truncate A p to have two components only, w 2 and 

W4 for instance, then p is the product of e~ r and a polynomial V(r), 

6 7 P{r) = 8r 6 - 12r 4 + 18r 2 + 9 + A p 2 12 (2t r)* 3“5 (2r 2 - 1) + Ap A 12 (tt/ 2)* 15“* (8r 4 - 14r 2 + 1). (20) 

Rescale out inessential factors, for a simpler polynomial, P = 8 r 6 — 12r 4 + 18r 2 + 9-|- Ai? 2 (2r 2 —1) + Ai? 4 ( 8 r 4 — 14r 2 + l). 
Eliminate r between P and dP/dr. The resultant 1Z(AR 2 , AR 4 ), when it vanishes, gives the border of the convex 
domain of parameters AR 2 , AR 4 where P remains positive definite. This domain contains the origin, because of 
p(r, 0). The precise form of 1Z is a little cumbersome and does not need to be published here. But the resulting border 
is shown in Figure 10. Generalizations to more A p parameters are obvious, with more cumbersome resultants TZ. 



FIG. 10. Domain of values of AR 2 and AR 4 acceptable for the positivity of the density of the toy model. The domain sits 
inside the full line curve and left of the straight line. It contains the origin. 
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VII. DISCUSSION AND CONCLUSION 


The subject of orthogonal polynomials has been so treated and overtreated that any claim to novelty must contain 
much more than a change of the integration measure. We took therefore a different approach, motivated by a law of 
physics and/or chemistry, matter conservation. This means a constraint of a vanishing average for the states described 
by weighted polynomials. 

For a support [0, oo[ and a simple exponential weight such as e ~^ r , a non trivial generalization of Laguerre poly¬ 
nomials occurs. This extends the generalization of Hermite polynomials described in [1] with the support ] — oo, oo[ 
and Gaussian weights such as e~i r . 

We also took care of cylindrical and spherical geometries, by replacing J dr with f dr r and J drr 1 2 3 4 5 6 7 , respectively. 
The new sets of constrained polynomials are clearly sensitive to the geometry. 

For finite supports such as [0,1] and constant weights, the constraint is already satisfied by the usual brand of 
orthogonal polynomials as soon as their order is > 1. In that sense, we did not find significantly original generalizations 
of Legendre polynomials, although we generated polynomials fitted to the cylindrical and spherical geometries. The 
cause of the failure is transparent: when the weight /i(r) is a constant, there is no difference between the orthogonality 
metric fi 2 and the constraint weight fi. 

For each set of new polynomials we found a recursion relation and a differential equation. There seems to be a 
systematic property for those cases where the constraint generates truly original polynomials, namely when /.t 2 ^ fi. 
In such cases, recursion and differentiation seem to be necessarily entangled. This does not happen for traditional 
orthogonal polynomials, indeed, and this “entanglement” may deserve some future attention. 

Constrained polynomials expressing matter conservation in centrifuges do make an original set if the fluid under 
centrifugation is compressible; a non constant reference weight /.i is indeed in order there. But the set depends on 
the precise form of /a via potentially many physical parameters. We found it difficult to design, through scaling, a 
sufficiently “universal” set. “Centrifuge polynomials” will have to be calculated specifically for each practical situation. 

For those new polynomials generalizing the Hermite and Laguerre ones, we found a description of the subspace 
accounting for their defect of completeness. A codimension 1 is the consequence of the constraint, expressed at first 
by the obvious lack of a polynomial of order n = 0. 

Finally the use of such polynomials was illustrated by a toy model for the Hohenberg-Kohn functional. A slightly 
surprising result was found: our polynomials, those of low order at least, define potential perturbations which are 
reflected by density perturbations having almost the same shapes. This occurs despite the delocalization created by 
the kinetic energy operator, hints at short ranges in effective interactions and validates the localization spirit of the 
Thomas-Fermi method. Whether such hints are good when the full zoology of the density functional is investigated 
is, obviously, an open question; for a review of the richness of the functional, we refer to [7]. If long range forces are 
active, a significant amount of delocalization between the “potential cause” and the “density effect” is not excluded. It 
would be interesting indeed to discover collective degrees of freedom in this connection between potential and density. 
In any case, our main conclusion may be that the new polynomials provide, for the context of matter conservation, a 
discrete and full set of modes and coordinates, hence a systematic and constructive representation of phenomena. 

It is a pleasure to thank Y. Abe, J.-P. Boujot, B. Eynard, C. Normand, R. Peschanski and A. Weiguny for stimulating 
discussions. 


[1] B.G. Giraud, M.L. Mehta and A. Weiguny, C. R. Physique 5 (2004) 781 

[2] P. Hohenberg and W. Kohn, Phys. Rev. 136 3B (1964) B864 

[3] N. D. Mermin, Phys. Rev. 137 5A (1965) A1441 

[4] R. Berg and L. Wilets, Proc. Phys. Soc. A, LXVIII, (1955) 229; W. Kohn and L. J. Sham, Phys. Rev. 140 4A (1965) A1133 

[5] A. Messiah, Mecanique Quantique, Dunod ed. Paris (1962) p. 362 

[6] J.-P. Boujot, These de 3e cycle, Calcul de moments appliques a 1’evaluation numerique d’integrales, Faculte des sciences de 
l’universie de Paris (1964) 

[7] J.P. Perdew and S. Kurth, in A Primer in Density Functional Theory, C. Fiolliais, F. Nogueira and M. Marques, eds., 
Lecture Notes in Physics, 620, Springer, Berlin (2003) 


11 




